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Abstract. We define the parametric closure problem, in which the input 
is a partially ordered set whose elements have linearly varying weights 
and the goal Is to compute the sequence of minimum-weight lower sets of 
the partial order as the weights vary. We give polynomial time solutions 
to many important special cases of this problem including semiorders, 
reachability orders of bounded-treewidth graphs, partial orders of bounded 
width, and series-parallel partial orders. Our result for series-parallel 
orders provides a significant generalization of a previous result of Carlson 
and Eppstein on bicriterion subtree problems. 


1 Introduction 


Parametric optimization problems are a variation on classical combinatorial 
optimization problems such as shortest paths or minimum spanning trees, in 
which the input weights are not fixed numbers, but vary as functions of a 
parameter. Different parameter settings will give different weights and different 
optimal solutions; the goal is to list these solutions and the intervals of parameter 
values within which they are optimal. As a simple example, consider maintaining 
the minimum of n input values as a parameter controlling these values changes. 
This parametric minimum problem asks for the lower envelope of a collection of 
input functions; for linear functions this is equivalent by projective duality to a 


planar convex hull 25 , and can be constructed in time 0(n log n); more general 
function classes such as piecewise-polynomial functions also have efficient lower- 
envelope algorithms [16| . The parametric minimum spanning tree problem (with 
linear edge weights) has polynomially many solutions that can be constructed in 
polynomial time ; the parametric shortest path problem has a number of 

solutions and running time that are exponential in log^ n on n-vertex graphs . 

As well as the obvious applications of this formulation to real-world problems 
with time-varying but predictable data (such as rush-hour route planning), 
parametric optimization problems have another class of applications, to bicriterion 
optimization. In bicriterion problems, each input has two numbers associated with 
it, that can be summed over the elements of a candidate solution. For instance, 
these two numbers might be the x and y coordinates of points in the plane, 
the mean and variance of a normal distribution, an initial investment cost and 
expected profit of a business opportunity, or the cost and log-likelihood of failure 
of a communications link. The goal is to find a solution that optimizes a nonlinear 
combination of these two sums of values, such as the distance from the origin, 
probability of exceeding a given threshold, percentage return on investment (the 
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Fig. 1. An instance of the parametric closure problem. Left: The Hasse diagram of a 
partially ordered set N of four elements, each with a weight that varies linearly with a 
parameter p. Center: The distributive lattice lower(A’) of lower sets of N. Right: The 
point set project(lower(A^)) and its convex hull. The upper hull (dashed) gives in left-to- 
right order the sequence of six distinct maximum-weight closures as the parameter p 
varies continuously from —oo to -|-(xj. 


ratio of total profit to total investment cost), or cost-reliability ratio. Many 
natural bicriterion problems can be expressed as finding the maximum of a 
quasiconvex function of the two sums (a function whose lower level sets are 
convex sets) or equivalently as finding the minimum of a quasiconcave function of 
the two sums. When this is the case, the optimal solution can always be obtained 
as one of the solutions of a parametric problem, defined by re-interpreting the 
two numbers associated with each input element as the slope and y-intercept 
of a linear parametric function j6,21 . In this way, any algorithm for solving a 
parametric optimization problem can also be used to solve bicriterion versions of 
the same type of optimization problem. 

In this paper we formulate and provide the first study of the parametric closure 
problem, the natural parametric variant of a classical optimization problem, the 
closure problem [^,28 . A closure in a directed graph is a subset of 


maximum 


vertices such that all edges from a vertex in the subset go to another vertex in the 
subset; the maximum closure problem is the problem of hnding the highest-weight 
closure in a vertex-weighted graph. Equivalently we seek the highest-weight lower 
set of a weighted partial order, where a lower set is a subset of the elements 
of a partial order such that if a: < y in the order, and y belongs to the subset, 
then X also belongs to the subset. Applications of this problem include open 
pit mining |23|, military attack planning |27|, freight depot placement |3 30 


scheduling with precedence constraints 8 22 , image segmentation j2 15 , stable 
marriage with maximum satisfaction 19 , and treemap construction in information 
visualization j^. Maximum closures can be found in polynomial time by a 
reduction to maximum flow [T7|[^ or by direct algorithms . 

In the parametric closure problem, we assign weights to the vertices of a 
directed graph (or the elements of a partial order) that vary linearly as functions 
of a parameter, and we seek the closures (or lower sets) that have maximum 
weight as the parameter varies. As described above, an algorithm for this problem 
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can also solve bicriterion closure problems of maximizing a quasiconvex function 
(or minimizing a quasiconcave function) of two sums of values. Although we have 
not been able to resolve the complexity of this problem in the general case, we 
prove near-linear or polynomial complexity for several important special cases of 
the parametric closure problem. 

We do not know of previous work on the general parametric closure prob¬ 
lem, but two previous papers can be seen in retrospect as solving special cases. 
Lawler 22 studied scheduling to minimize weighted completion time with prece¬ 


dence constraints. He used the closure that maximizes the ratio x/y of the priority 
X and processing time ?/ of a job or set of jobs to decompose instances of this 
problem into smaller subproblems. As Lawler showed, the optimal closure can 
be found in polynomial time by a binary search where each step involves the 
solution of a weighted closure problem. Replacing the binary search by parametric 
search 26 would make this algorithm strongly polynomial; however, both search 
methods depend on the specific properties of the ratio function and would not 
work for other bicriterion problems. A second paper, by Carlson and Eppstein [^, 
considers bicriterion versions of the problem of finding the best subtree (contain¬ 
ing the root) of a given rooted tree with weighted edges. These subtrees can be 
modeled as lower sets for a partial order on the tree edges in which two edges are 
comparable when they both belong to a path from the root; this partial order is 
series-parallel, and we greatly generalize the results of Carlson and Eppstein in 
our new results on parametric closures for series-parallel partial orders. 


Parametric optimization as an implicit convex hull problem. Parametric 
optimization problems can be formulated dually, as problems of computing convex 
hulls of implicitly defined two-dimensional point sets. Suppose we are given a 
parametric optimization problem in which weight of element i is a linear function 
OiA -f of a parameter A, and in which the weight of a candidate solution S (a 
subset of elements, constrained by the specific optimization problem in question) 
is the sum of these functions. Then the solution value is also a linear function, 
whose coefficients are the sums of the element coefficients: 

A+ 

iGS / ViGS 

Instead of interpreting the numbers and hi as coefficients of linear functions, we 
may re-interpret the same two numbers as the x and y coordinates (respectively) 
of points in the Euclidean plane. In this way any family T of candidate solutions 
determines a planar point set, in which each set in T corresponds to the point 
given by the sum of its elements’ coefficients. We call this point set project(J^), 
because the sets in T can be thought of as vertices of a hypercube Qn = {0,1}" 
whose dimension is the number of input elements, and project determines a linear 
projection from these vertices to the Euclidean plane. 

Let hull(project(J^)) denote the convex hull of this projected planar point 
set. Then for each parameter value the set in J- minimizing or maximizing the 
parameterized weight corresponds by projective duality to a vertex of the hull. 


^ ' OjA -j- hi — 
jGS 
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and the same is true for the maximizer of any quasiconvex function of the two 
sums of coefficients and bi- Thus, parametric optimization can be reformulated 
as the problem of constructing this convex hull, and bicriterion optimization can 
be solved by choosing the best hull vertex. 

New results. For an arbitrary partially ordered set P, define lower(P) to be the 
family of lower sets of P. As a convenient abbreviation, we define polygon(P) = 
hull(project(lower(P))). We consider the following classes of partially ordered set. 
For each partial order P in one of these classes, we prove polynomial bounds on 
the complexity of polygon(P) and on the time for constructing the hull. These 
results imply the same time bounds for parametric optimization over P and for 
maximizing a quasiconvex function over P. 

Semiorders. This class of partial orders was introduced to model human pref¬ 
erences |24| in which each element can be associated with a numerical value, 
pairs of elements whose values are within a fixed margin of error are incompa¬ 
rable, and farther-apart pairs are ordered by their numerical values. For such 
orderings, we give a bound of 0(n logn) on the complexity of polygon(P) and 
we show that it can be constructed in time 0(nlog^ n) using an algorithm 
based on the quadtree data structure. 

Series-parallel partial orders. These are orders formed recursively from smaller 
orders of the same type by two operations: series compositions (in which all 
elements from one order are placed earlier in the combined ordering than 
all elements of the other order) and parallel compositions (in which pairs 
of one element from each ordering are incomparable). These orderings have 
been applied for instance in scheduling applications by Lawler |22| . For such 
orderings, the sets of the form polygon(P) have a corresponding recursive 
construction by two operations: the convex hull of the union of two convex 
polygons, and the Minkowski sum of two convex polygons. It follows that 
polygon(P) has complexity 0{n). This construction does not immediately 
lead to a fast construction algorithm, but we adapt a splay tree data structure 
to construct polygon(P) in time O(nlogn). Our previous results for optimal 
subtrees follow as a special case of this result. 

Bounded treewidth. Suppose that partial order P has n elements and its 
transitive reduction forms a directed acyclic graph whose underlying undi¬ 
rected graph has treewidth w. (For prior work on treewidth of partial orders, 
see |20].) Then we show that polygon(P) has polynomially many vertices, 
with exponent 0{w), and that it can be constructed in polynomial time. 
Incidence posets The incidence poset of a graph G has the vertices and edges 
as elements, with an order relation x < y whenever x is an endpoint of y. 
One of the initial applications for the closure problem concerned the design 
of freight delivery systems in which a certain profit could be expected from 
each of a set of point-to-point routes in the system, but at the cost of setting 
up depots at each endpoint of the routes j^[^; this can be modeled with 
an incidence poset for a graph with a vertex at each depot location and 
an edge for each potential route. Since the profits and costs have different 
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timeframes, it is reasonable to combine them in a nonlinear way, giving a 
bicriterion closnre problem. The transitive reduction of an incidence poset 
is a subdivision of G with the same treewidth, so our technique for partial 
orders of bounded treewidth also applies to incidence posets of graphs of 
bounded treewidth. 

Bounded width. The width of a partial order is the maximum number of 
elements in an antichain, a set of mutually-incomparable elements. Low- 
width partial orders arise, for instance, in the edit histories of version control 
repositories |^. The treewidth of a partial order is at most equal to its width, 
but partial orders of width w have Oin^) lower sets, tighter than the bound 
that would be obtained by using treewidth. We show more strongly using 
quadtrees that in this case polygon(P) has -l-nlogn) vertices and 

can be constructed in time within a logarithmic factor of this bound. 

We have been unable to obtain an example of a family of partial orders with 
a nonlinear lower bound on the complexity of polygon(P), nor have we been able 
to obtain a nontrivial upper bound on the hull complexity for unrestricted partial 
orders. Additionally, we have been unable to obtain polynomial bounds on the 
hull complexity of the above types of partial orders for dimensions higher than 
two. We also do not know of any computational complexity bounds (such as 
NP-hardness) for the parametric closure problem for any class of partial orders 
in any finite dimension. We leave these problems open for future research. 

For space reasons we describe only the semiorder and series-parallel results in 
the main text of our paper, deferring the remaining results to appendices. 

2 Minkowski sums and hulls of unions 

Our results on the complexity of the convex polygons polygon(P) associated with 
a partial order hinge on decomposing these polygons recursively into combinations 
of simpler polygons. To do this, we use two natural geometric operations that 
combine pairs of convex polygons to produce more complex convex polygons. 

Definition 1. For any two convex polygons P and Q, let P (B Q denote the 
Minkowski sum of P and Q (the set of points that are the vector sum of a point 
in P and a point in Q), and let P ^ Q denote the convex hull of the union of P 
and Q. 

Lemma 1 (folklore) . If convex polygons P and Q have p and q vertices respec¬ 
tively, then P 0 (5 and P dJ Q have at most p-\- q vertices, and can he constructed 
from P and Q in time 0(p + q). 

We omit the (easy) proof for space reasons. 

Corollary 1. Suppose that P is a convex polygon, described as a formula that 
eombines a set of n points in the plane into a single polygon using a sequenee 
of 0 and y operations. Suppose in addition that, when written as an expression 
tree, this formula has height h. Then P has at most n vertiees and it may be 
constructed from the formula in time 0{nh). 
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More complex data structures can reduce this time to 0{n log n); see Section]^ 
In higher dimensions, the convex hull of n points and Minkowski sum of n 
line segments both have polynomial complexity with an exponent that depends 
linearly on the dimension. However, we do not know of an analogous bound on the 
complexity of convex sets formed by combining Minkowski sum and hull-of-union 
operations. If such a bound held, we could extend our results on parametric 
closures to the corresponding higher dimensional problems. 

3 Semiorders 

A semiorder is a type of partial order defined by Luce [24] to model human 
preferences. Each element of the order has an associated numerical value (its 
utility to the person whose preferences are being modeled). For items whose 
utilities are sufficiently far from each other, the ordering of the two items in 
the semiorder is the same as the numerical ordering of their utilities. However, 
items whose utilities are within some (global) margin of error of each other are 
incomparable in the semiorder. Similar concepts of comparisons of numerical 
values with margins of error give rise to semiorders in many other areas of science 
and statistics [^. For efficient computations on semiorders we will assume that 
the utility values of each element are part of the input to an algorithm, and that 
the margin of error has been normalized to one. For instance, the semiorder N of 


b, c, and d respectively. With this information in hand, the comparison between 
any two elements can be determined in constant time. 

The concept of a lower set is particularly natural for a semiorder: it is a set of 
elements whose utility values could lie below a sharp numerical threshold, after 
perturbing each utility value by at most half the margin of error. In this way, 
the closure problem (the problem of finding a maximum weight lower set) can 
alternatively be interpreted as the problem of finding the maximum possible 
discrepancy of a one-dimensional weighted point set in which the location of each 
point is known imprecisely. Semiorders may have exponentially many lower sets; 
for instance, if all items have utilities that are within one unit of each other, all 
sets are lower sets. Nevertheless, as we show in this section, if S' is a semiorder, 
then the complexity of polygon(S) is near-linear. 

If S is any parametrically weighted semi-order, we may write the sorted order 
of the utility values of elements of S as uq,ui, ..., u„_i where n = |S|, and the 
elements themselves (in the same order) as xq, xi,..., Xn-i- By padding S with 
items that have a fixed zero weight and a utility that is smaller than that of 
the elements by more than the margin of error, we may assume without loss of 
generality that n is a power of two without changing the values of the parametric 
closure problem on S. 

Definition 2. Let L be an arbitrary lower set in lower(S). Let j be the largest 
index of an element Xi of L. Let i be the smallest index of an element Xi sueh 
that Xi does not belong to L and i < j, or —1 if no such element exists. Define 
extremes(L) to he the pair of integers {i + l,j). 


Figure 1 can be represented as a semiorder with utilities 2/3, 0, 2, and 4/3 for a, 
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Fig. 2. The grid [0,n — 1]^, with the two regions that cannot be part of the image 
of extremes. The left image shows a square subproblem s and free(s); the right image 
shows the quadtree decomposition of the grid used to prove [Theorem 1[ 


Thus, extremes maps the family lower(S') to the integer grid [0, n — 1]^, with 
potentially many lower sets mapped to each grid point. However, not every grid 
point is in the image of lower(5'): a point {i,j) with i > j cannot be the image of 
a lower set, because the element defining the first coordinate of extremes must 
have an index smaller than the element defining the second coordinate. And when 
i > 0, a point (i,j) with Ui-i < Uj — 1 (i.e. with utility values that are beyond 
the margin of error for the semiorder) also cannot be the image of a lower set, 
because in this case Xi-i < Xj in the semiorder, so every lower set that includes 
Xj also includes Xi-i. Thus, the image of extremes lies in an orthogonally convex 
subset of the grid, bounded below by its main diagonal and above by a monotone 
curve (Figure 2). 


Definition 3. Let s be any square subset of the integer grid [0, n— 1]^, and define 
subproblem(s) to be the partially-ordered subset of the semiorder S consisting of 
the elements whose indices are among the rows and columns of s. Define free(s) 
to be the (unordered) set of elements of S that do not belong to subproblem(s), 
but whose indices are between pairs of indices that belong to subproblem(s). (See 
Figure left, for an example.) 


Observation 1 Given a square s, suppose that the subfamily T of lower(S') that 
is mapped by extremes to s is nonempty. Then each set in T is the disjoint union 
of a lower set o/subproblem(s) and an arbitrary subset o/free(s), and all such 
disjoint unions belong to J-. 


Observation 2 For any square s, let powerset(free(s)) be the family of all subsets 
o/free(s). Let weight function w : S define a projection project from families 

of sets to point sets in Then project(powerset(free(s))) is the Minkowski 
sum of the sets {(0,0),r(;(a;i)} for Xi S free(s). Its convex hull is a centrally 
symmetric convex polygon hull(proJect(powerset(free(s)))) (the Minkowski sum of 
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the corresponding line segments) with at most k = 2|free(s)| sides, and can he 
constructed in time 0{k\ogk). 

Corollary 2. Given a square s, letiF be the sub-family o/lower(iS') that is mapped 
by extremes into s, and assume a weight function w defining a projection project. 
Then hull(proJect(J^)) = polygon(subproblem(s)) 0 hull(project(powerset(free(s)))). 

Lemma 2. Let s be a square in the grid [0,n — 1]^, and subdivide s into four 
congruent smaller squares Si (0 < i < 4). Let polygon (subproblem(s)) have c 
vertices, define Ci in the same way for each Si, and let t be the side length of s. 
Then c < ^ Ci + 0{£), and polygon(subproblem(s)) can be constructed from the 
corresponding hulls for the smaller squares in time Oif^ Ci + £\ogt). 

Proof. For each smaller square Si, define Ti to be the subset of subproblem(s) 
mapped to Si, and define Hi = hull(proJect(J'i)). Then by Corollary(viewing 
Si as a subproblem of subproblem(s)), 

Hi = polygon(subproblem(s)) 0 hull(proJect(powerset(free(s) \ free(si))). 

The set free(s) \ free(si) has cardinality 0[t), so by Lemmaand Observation 
Hi has complexity Ci 0 0{£) and can be constructed in time 0{ci 0 £\ogt). 
Applying Lemma again, polygon(subproblem(s)) = Hq dJ i/i dJ H 2 H 2 has 
complexity at most X) + 0{£) and can be constructed from the polygons Hi in 
time 0(^ Ci 0 £ log £). □ 

Theorem 1. If S is a semiorder with n elements Xi, specified with their utility 
values Ui and a system of two-dimensional weights w{xi), then polygon(S') has 
complexity O(nlogn) and can he constructed in time O(nlog^n). 

Proof. We sort the utility values, pad n to the next larger power of two if necessary 
and form a quadtree decomposition of the grid [0,n — 1]^ (as shown in 
right). For each square s of this quadtree, we associate a convex polygon (or 
empty set) polygon(subproblem(s)) computed according to the following cases: 

— If s is a subset of the grid points for which i > j, or for which Ui-i < Uj — 1, 
then no lower sets are mapped into s by extremes. We associate square s with 
the empty set. 

— If s is a subset of the grid points for which i < j and Ui-i > Uj — 1, then every 
two elements of subproblem(s) are incomparable. In this case, we associate 
square s with the polygon hull(proJect(powerset(subproblem(s)))) computed 
according to Observation 

— Otherwise, we split s into four smaller squares. We construct the polygon 
associated with s by using Lemma to combine the polygons associated with 
its children. 

It follows by induction that the total complexity of the polygon constructed at 
any square s of the quadtree is OC^ii), and the total time for constructing it 
is 0(^£ilog£i)) where £i is the side length of the ith square of the quadtree 


Figure 2 



and the sum ranges over all descendants of s. As a base case for the induction, 
a square containing only a single grid point is associated with a subproblem 
with one element, with only one lower set that maps to that grid point, and a 
degenerate convex polygon with a single vertex. The polygon constructed at the 
root of the quadtree is the desired output, and it follows that it has combinatorial 
complexity and time complexity of the same form, with a sum ranging over all 
quadtree squares. 

The conditions i > j and Ui-i < Uj — 1 define two monotone curves through 
the grid, and we split a quadtree square only when it is crossed by one of these 
two curves. It follows that the squares of side length £ that are subdivided as part 
of this algorithm themselves form two monotone chains, and that the number 
of all squares of side length £ is 0(n/£). The results of the theorem follow by 
summing up the contributions to the polygon complexity and time complexity 
for the O(logn) different possible values of £. □ 


4 Series-parallel partial orders 


Series-parallel partial orders were considered in the context of a scheduling 
problem by Lawler 22 , and include as a special case the tree orderings previously 


studied in our work on bicriterion optimization [^. They are the partial orders 
that can be constructed from single-element partial orders by repeatedly applying 
the following two operations: 


Series composition. Given two series-parallel partial orders Pi and P 2 -, form 
an order from their disjoint union in which every element of Pi is less than 
every element of P 2 . 

Parallel composition. Given two series-parallel partial orders Pi and P 2 , form 
an order from their disjoint union in which there are no order relations 
between Pi and P 2 . 


Observation 3 If P is the series composition of Pi and P 2 , then polygon(P) is 
the convex hull of the union of polygon (Pi) and a translate (by the sum of the 
weights of the elements of Pi) 0 / polygon (P 2 ). If P is the parallel composition of Pi 
and P 2 , then polygon(P) is the Minkowski sum 0 / polygon(Pi) and polygon(P 2 ). 

Recursively continuing this decomposition gives us a formula for polygon (P) 
in terms of the itu and oplus operations. By Lemmaj^we immediately obtain: 

Corollary 3. If P is a series-parallel partial order with n elements, then polygon(P) 
has at most 2n vertices. 


However, the depth of the formula for polygon(P) may be linear, so using 
Lemma to construct polygon(P) could be inefficient. We now describe a faster 
algorithm. The key idea is to follow the same formula to build polygon(P), but 
to represent each intermediate result (a convex polygon) by a data structure 
that allows the ItU and 0 operations to be performed more quickly for pairs of 
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polygons of unbalanced sizes. Note that a Minkowski sum operation between a 
polygon of high complexity and a polygon of bounded complexity can change a 
constant fraction of the vertex coordinates, so to allow fast Minkowski sums our 
representation cannot store these coordinates explicitly. 


Lemma 3. It is possible to store convex polygons in a data structure such that 
destructively merging the representations of two polygons of m and n vertices 
respectively by a ^ or 0 operation (with m < n) can be performed in time 
0(mlog((m + n)/m)). 

Proof. We store the lower and upper hulls separately in a binary search tree data 
structure, in which each node represents a vertex of the polygon, and the inorder 
traversal of the tree gives the left-to-right order of the vertices. The node at the 
root of the tree stores the Cartesian coordinates of its vertex; each non-root node 
stores the vector difference between its coordinates and its parents’ coordinates. 
Additionally, each node stores the vector difference to its clockwise neighbor 
around the polygon boundary. In this way, we can traverse any path in this tree 
and (by adding the stored vector difference) determine the coordinates of any 
vertex encountered along the path. We may also perform a rotation in the tree, 
and update the stored vector differences, in constant time per rotation. 

We will keep this tree balanced (in an amortized sense) by using the splay 
tree balancing strategy |31] : whenever we follow a search path in the tree, we 
will immediately perform a splay operation that through a sequence of double 
rotations moves the endpoint of the path to the root of the tree. By the dynamic 
finger property for splay trees , a sequence of m accesses in sequential order 
into a splay tree of size n will take time 0{m\og{(rn + n)/m)). 

To compute the hull of the union (the ity operation) we insert each vertex 
of the smaller polygon (by number of vertices), in left-to-right order, into the 
larger polygon. To insert a vertex v, we search the larger polygon to find the 
edges with the same x-coordinate as v and use these edges to check whether v 
belongs to the lower hull, the upper hull, or neither. If it belongs to one of the 
two hulls, we search the larger polygon again to find its two neighbors on the 
hull. By performing a splay so that these neighbors are rotated to the root of the 
binary tree, and then cutting the tree at these points, we may remove the vertices 
between v and its new neighbors from the tree without having to consider those 
vertices one-by-one. We then create a new node for v and add its two neighbors 
as the left and right child. 

To compute the Minkowski sum (the © operation) we must simply merge 
the two sequences of edges of the two polygons by their slopes. We search for 
each edge slope in the smaller polygon. When its position is found, we splay the 
vertex node at the split position to the root of the tree, and then split the tree 
into its left and right subtrees, each with a copy of the root node. We translate 
all vertices on one side of the split by the vector difference for the inserted edge 
(by adding that vector only to the root of its tree), and rejoin the trees. □ 
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Theorem 2. If P is a series-parallel partial order, represented by its series- 
parallel deeomposition tree, then polygon (P) has eomplexity 0(n) and may he 
construeted in time O(nlogn). 

Proof. We follow the formula for constructing polygon (P) by ItU and 0 operations, 
using the data structure of Lemma We charge each merge operation to the 
partial order elements on the smaller side of each merge. If a partial order 
element belongs to subproblems of sizes uq = 1, ni, Uh = n where h is 
the height of the element, then the time charged to it is 0{\og{ni/ni-i)) = 
0(logn*(«*/’^i-i)) = O(logn). □ 

Acknowledgements. Supported in part by NSF grant 1228639 and ONR grant 
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A Bounded tree-width 


In this section, we consider partial orders defined by reachability on directed 
acyclic graphs whose underlying undirected graph has bounded tree width. 

A tree-decomposition of an undirected graph G is a tree T whose vertices 
are called bags, each of which is associated with a set of vertices of G, with the 
following two properties: 

— Each vertex of G belongs to a set of bags that induces a connected subtree 

of T, and 

— For each edge of G there exists a bag containing both endpoints of the edge. 

The width of a tree-decomposition is one less than the maximum number of vertices 
in a bag, and the treewidth of G is the minimum width of a tree-decomposition 
of G. The treewidth, and an optimal tree-decomposition, can be found in time 
linear in the number of vertices of a given graph and exponential in its width. 

By splitting nodes with high degree into more than one node (all associated 
with equal sets of vertices of G), we may assume without loss of generality that 
T has maximum degree three. We may also assume without loss of generality 
that T has 0{n) nodes. These assumptions do not change the width of the 
decomposition. 

We will use the following well-known result: 

Lemma 4 (folklore). Any k-node forest F of maximum degree three can he 
partitioned by the deletion of one node x into two smaller forests, each having at 
most 2k/i nodes. 

Proof. Let x be chosen to minimize the maximum size of a tree in the remaining 
forest. Then each tree has at most 2fc/3 nodes: otherwise we could reduce the 
maximum tree size by moving x into the largest tree or (if it is already there) 
moving x to one of its neighbors. 

We use a greedy strategy to select trees for one of the two forests, largest 
first, until the forest has at least k/Z nodes. If this happens on the first step, we 
are done; otherwise, each step increases the size of the forest by less than fc/3 
nodes, so the size of the forest cannot jump from less than k/2> to greater than 
2k/?> at any step. □ 

Lemma 5. Let P be a partial order with n elements whose underlying graph has 
treewidth w. Then polygon (P) can be represented as a formula of binary lyj and 
0 operations with height at most (w + 2 )(log 3^2 

Proof. We recursively decompose the given tree decomposition. At each step of 
the decomposition we will have a subforest F of the original tree decomposition 
(initially T itself), a bag B of F chosen according to Lemma (initially, the 
choice given by the lemma for T), and an assignment of some elements of P 
to two disjoint subsets L and U of P (initially empty). The meaning of this 
information is that we restrict our attention to lower sets on the elements of P 
that participate in bags of F, and that are consistent with including the elements 
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of L in the lower set and excluding the elements of U from the lower set. Each 
subproblem is associated with a polygon, the convex hull of the projections of 
the lower sets it represents. If an element of L is greater than an element of U, 
then there are no consistent lower sets and the associated polygon is empty. 

Given this information, suppose that some element b oi B does not already 
belong to LU U. In this case, we partition the problem into two subproblems: 
one in which b has been added to L and the other in which b has been added 
to U. The polygon representing the problem prior to this choice of which set to 
place b into can be represented as the convex hull of the union of the polygons 
for these two subproblems, because they represent disjoint families of lower sets 
that together cover all the lower sets represented by the whole problem. 

Alternatively, suppose that the elements of bag B are completely covered by 
LU U. In this case, let Fi and F 2 be the two subforests of F, of size smaller 
by a factor of 2/3 or better, guaranteed to exist by Lemma We can form 
subproblems for Fi and F 2 by using Lemma to choose a bag within each of 
these two forests. Because all of the elements of bag B have already been placed 
into L or into U, all remaining choices about the membership of the lower sets in 
Fi and in F 2 are completely independent of each other. Therefore, the polygon 
representing the given problem can be represented as the Minkowski sum of the 
polygons representing Fi and L 2 . 

There can be at most log 3/2 n + 0(1) steps of this recursive decomposition at 
which we split the current forest F into two smaller subforests, because each such 
step reduces the forest size by a 2/3 or better factor. Between any two such steps 
there can be at most w + 1 steps at which we choose whether to place a vertex of 
the currently chosen bag into L or U, because there are at most w + 1 vertices to 
place (fewer if the bag is smaller than the maximum bag size or shares vertices 
with a previously chosen bag). Therefore, the total height of the decomposition 
is as stated. □ 

By choosing the initial bag for each subproblem more carefully, we can also 
represent polygon(P) as an acyclic circuit of dJ and 0 operations with 0(2’"n) 
complexity and 0(z(;logn) depth (but with a somewhat higher constant factor 
in the depth). However, we omit the details as we have been unable to generalize 
Corollary [2 from formulas to circuits. 

Theorem 3. Let P be a partial order with n elements whose underlying graph 
has treewidth w. Then polygon (P) has eomplexity and can be constructed 

in time . 

Proof. By Lemmaj^ there exists a formula of height {w + 2)(log3/2 n + 0(1)) for 
polygon (P). Because this is a formula of binary operations, the total size of the 
formula is 

2(u)+2)(log3y2n+0(l)) _ 

The result follows by applying Corollary to this formula. □ 

Corollary 4. Let P be the incidence poset of an n-vertex graph of treewidth w. 
Then polygon(P) has complexity and can be constructed in time 
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For generalized fences and polytrees, we obtain a polynomial bound on the 
complexity of polygon(P) by applying Theoremwith w = 1, but we can tighten 
this bound in three ways. First, by applying the recursive decomposition to the 
given tree rather than to its tree decomposition, we need only choose whether 
to place into L oi U a single element for each tree node, rather than the two 
elements of a bag. Second, we observe that there are only 0{n) sequences of 
choices of which of two forests to continue the recursion in, because each such 
sequence ends at a single tree node; thus, we can charge all of these choices to a 
single 0(n) factor in the complexity rather than charging a factor of two for each 
level of formula depth given by these choices. And finally, when the underlying 
polytree is a path (as it is in the fences and in some definitions of generalized 
fences) we may always split each subproblem at the path midpoint, reducing the 
log 3/2 terms in the previous analysis to log 2 n. Based on these observations, we 
obtain: 

Corollary 5. Let P be the reachability order of a polytree T. Then polygon(P) 
has complexity and can be constructed in time logn). 

IfT is an oriented path, polygon(P) has complexity 0{n^) and can be constructed 
in time O(nlogn). 
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B Bounded width 


The width of a partially ordered set is the size of its largest independent set. 
Partially ordered sets of bounded width arise, for instance, in the version histories 
of a distributed version control repository controlled by a small set of developers; in 
this application, there may be many elements of the partially ordered set (versions 
of the repository) but the width is bounded by the number of developers . A 
lower set in this application is a set of versions that could possibly describe the 
simultaneous states of all the developers at some past moment in the history of 
the repository. 

The lower sets of a partially ordered set correspond one-for-one with its 
independent sets: each lower set is uniquely determined by an independent set, 
the set of maximal elements of the lower set. Therefore, in a partially ordered set 
of width w with n elements there can be at most 0{n'^) lower sets. More precisely, 
by Dilworth’s theorem every partial order of width w can be partitioned 
into w totally-ordered subsets (chains); the number of lower sets can be at most 
the product of the numbers of lower sets of these chains, {n/w + 1)“'. As we show, 
however, there is a tighter bound for the complexity of polygon(P). 

We consider first the case w = 2. The more general case of bounded width 
will follow by similar reasoning. Thus, let P be a partially ordered set with n 
elements and width two. By Dilworth’s theorem, P can be decomposed into two 
chains, and every lower set L can be described by a pair of integer coordinates 
{x,y), where x is the number of elements of L that belong to the first chain and 
y is the number of elements of L in the second chain. 

Observation 4 For a width-two partial order P as described above, the set of 
pairs of coordinates {x, y) describing the lower sets of P forms an orthogonally 
convex subset of the integer grid [0,n]^, and the edges between points one unit 
apart in this set of grid points form the covering relation of the distributive lattice 
of lower sets. 

[Figure 3| depicts an example. We can use the following definition and observa¬ 
tion to partition the parametric closure problem for P into smaller subproblems 
of the same type. 

Definition 4. Let R be any grid rectangle. Then subproblem(P) is the family of 
lower sets of P whose grid points lie in R. 

Observation 5 For any grid rectangle, each set in subproblem(P) can be repre¬ 
sented uniquely as a disjoint union A\J B where A is the family of all elements 
of P whose coordinate value is below or to the left of R, and B is a lower set of 
the restriction of P to the elements whose coordinate value is within R. 

Observation 6 Suppose that every integer point of a grid rectangle R corre¬ 
sponds to a lower set of P. Then the restriction of P to the elements whose 
coordinate value is within R is a partial order in the form of the parallel compo¬ 
sition of two chains, a special case of a series-parallel partial order. R follows 
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Fig. 3. Hasse diagram of a width-two partial order partitioned into two chains (left) 
and the subset of the integer lattice formed by its lower sets (right) 



from Corollary^ that hull(project(subproblem(i?))) has complexity proportional 
to the perimeter of R. 

Theorem 4. If P is a partial order of width two, then polygon (P) has complex¬ 
ity 0{n\ogn) and can he constructed from the covering relation of P in time 
0{n\og^ n). 

Proof. From the covering relation, we can determine a partition into two chains 
and trace out the boundaries of the orthogonal grid polygon describing the lower 
sets of P, in linear time. We use a quadtree to partition the lower sets of P into 
squares, each with a side length that is a power of two, such that each grid point 
within each square corresponds to one of the lower sets of P. The squares of each 
size form two monotone chains within the grid, so the total perimeter of all of 
these squares is 0(n log n). 

For each of the power-of-two subintervals / of each of the two grid coordinates, 
we compute the convex hull of the lower sets of the corresponding chain whose 
grid points lie within /, as the hull of the union of the two previously-constructed 
polygons for the two halves of I. This computation takes time 0(n log n) over¬ 
all. After this step, each quadtree square’s polygon can be constructed as the 
Minkowski sum of the polygons for its two defining intervals, in time proportional 
to the side of the square. Therefore, the total time to compute the polygons for 
all of the quadtree squares is 0(n logn). 

It remains to form the convex hull of the union of these polygons. As the 
convex hull of O(nlogn) points, this takes time 0(n log^ n). □ 

For higher widths w, the same idea works (using an octree in three dimensions, 
etc). The total complexity of polygon(P) is proportional to the sum of side lengths 
of octree squares, 0(n’"“^), and the construction time is within a logarithmic 
factor of this bound. 
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